Quantum correlations in the collective spin systems 
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Quantum and classical pairwise correlations in two typical collective spin systems, i.e. the Dicke 
model and the Lipkin-Meshkov-Glick model are discussed. In the thermodynamical limit, analytical 
expressions of these correlations are derived. Within the effective technique proposed previously, 
theses correlations for finite size system can be also obtained numerically up to very large size. 
The scaling behavior for the quantum discord and its first derivative are analyzed. The comparison 
of the quantum discord and the concurrence are performed and essentially different behaviors are 
observed. 
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I. INTRODUCTION 

Correlation, as a fundamental feature, has been exten- 
sively investigated in many-body physics and the quan- 
tum information science [H, 0- They indeed unravel 
the key physical properties and characterize the remark- 
able phenomena of the critical systems, such as quantum 
phase transition (QPT) Q- QPTs reveal the qualitative 
change of the quantum systems resulted from the energy 
level crossings at zero temperature. 

Both the quantum and classical sources contribute to 
the correlation [J]. As one kind of the quantum correla- 
tion, the entanglement is widely applied as the character- 
istic trait to detect the quantum correlation Q. Specifi- 
cally, the entanglement successfully identifies the critical 
behaviors of the QPTs @-[I3< 

However, the entangle- 
ment may fail to capture the existence of the quantum 
correlation in some mixed separate states, in which the 
entanglement is considered not a good measure [H, EH ■ 
While a new kind of the quantum correlation, quantum 
discord (QD), provides the alternative route for measure- 
ment, which is present even under separable states [l2j]. 
The definition of the QD can be interpreted as the dif- 
ference of the total quantum information of the two sub- 
systems A and B before and after the local operation on 
the one of them. The QD has been proved as a good 
measure of the non-classical correlations beyond entan- 
glement. Furthermore, the QD has been indicated as the 
source to speed up the quantum computations [l4], [H| • 

Much attention has been recently paid to apply the QD 
to quantify the critic al prop erties of spin chain systems, 
especially the QPTs (lO-E^. Specifically, the QPTs of 
one dimensional XXZ model and transverse field Ising 
model have been anal ytic ally characterized by the QD at 
zero temperature (l6l. Il7l| and the finite temperature (T§l . 
Il9j . The critical properties of the half spin XY chain 



have been exhibited [20j. The QD has been also studied 
in extended XY models [2lil22|, chains. Apart from spin 
systems, the QD has been used to characterize the QPTs 
in the correlated electron systems [23| and a topological 
QPT in the Castelnovo-Chamon model[24j]. 



The Dicke model[25| and Lipkin-Meshkov-Glick 
(LMG) model[26[ are two well known quantum collective 
spin ones. The Dicke model describes the interaction of 
N two-level atoms with a single bosonic mode and ex- 
hibits QPT at the critical atom-cavity coupling. It can 
be realized experimentally in a Bose-Einstein condensate 
in an optical cavity 27]. LMG model was originally in- 
troduced in nuclear p hysics, but now has found appli- 
cations in other fields [28-30]. It describes N mutually 
interacting spins in a transverse magnetic field, and also 
undergos QPT and shares the same universality with the 
Dicke model. Except preliminarily approximate results 
for the QD in the LMG model pjj in the thermodynam- 
ical limit, the QD in these collective spin systems has not 
been well studied to date. 
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In the present paper we study the QD and the classi- 
cal correlation of the Dicke model and the LMG model 
both in the thermodynamic limit and for finite size sys- 
tems. The critical behavior related to the QD and its 
first derivative are analyzed. The paper is outlined as 
follows. In Sec. II we review the definition of the QD 
and the classical correlation, and then apply them to the 
spin collective models. In Sec. Ill we analyze the QD, 
the counterpart classical correlation, and the entangle- 
ment in detail of the Dicke model, and discussions are 
also presented. In Sec. IV, these features are studied in 
LMG. Finally, we summarize our work in Sec. V. 
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II. GENERAL FORMALISM FOR QUANTUM 
AND CLASSICAL CORRELATIONS IN THE 
COLLECTIVE SPIN MODELS 

The mutual information to describe the total correla- 
tion of the sub-systems A and B can be written as 



1(A; B) = H(A) + H(B) - H(A, B), 



(1) 



where H{a) = ~Y,kPa=k^ogp a =k (a = A,B) with 
p a =k the probability of the realization k for the sub- 
system a. The joint entropy of A and B is denoted as 

H(A,B) = -J2 3 ,kPA=],B=k\0gPA=j : B=k with PA=j,B=k 

the joint probability of the sub-systems A and B, re- 
alized by j and fc. By using the Bayes rule PA\B=k — 
PA,B=k/ps=k, the classical mutual information can also 
be rewritten into the equivalent expression 



J(A;B) = H(B)-H(A\B), 



(2) 



where H(A\B) = -J2^kPA= 3 \B=k^ogp A =j\B=k is the 
conditional entropy of the A and B. PA=j\B=k 1S the 
corresponding conditional probability. 

For the quantum systems, the two definitions of the 
mutual information in Eqs. (TT|) and are replaced 
by the Von Neumann entropy. The joint entropy is 
shown as H(A,B) = -^a.b{pa.B log Pa.b}, where p A ,B 
is the density matrix of the total system. Similarly, the 
Von Neumann entropy for the sub-system is H(A) = 
-Ti-AiPA^ogpA}, with pa = ^b{pa,b\- The quantum 
conditional entropy quantifies the missing information of 
A after selecting the state of B. It can be carried out 
from the conditional density matrix 



(3) 



where nf is the projector onto the state k of the sub- 
system B and the probability is pk — Tr^ e{nf Pa,b}- 
Usually, these two mutual information become different, 
which leads to the emergence of the quantum discord 
THA : B) = X(A;B) - J(A;B), described in Refs. [H, 
131( 1 . In the quantum measurements, the quantum discord 
between A and B reads 

V(A : B) = mm{H(A) H(A, B) + H(A\{U B k })}. (4) 
Then the classical correlation is described as 



C(A : B) = max{ff (A) 



H(A\{n B k })}. 



(5) 



Here, we apply the quantum discord to the collective 
spin systems, where the Dicke and LMG models are two 
typical examples. Since the quantum discord is adopted 
to describe the nonlocal quantum correlation, it is neces- 
sary to derive the pairwise density matrix. In this paper, 
we only study the pairwise correlations. For the Dicke 
model and LMG model, the reduced pairwise matrix in 



the standard basis ||' 
and a z \l) = -||)) 3: 



t),|tD,llt>,jU)}(witha z |t) = |D 
I , can be derived as 



P 



x + 
x + 

\u 



x + 

w 

V 



y 

w 

X- 



u, 

x* 
x* 

V- 



The detailed expressions for these elements are 

N 2 - 2N + A{Jl)±A{N - \){J Z ) 
AN(N - 1) ' 
(N-1)(J + )±([J + ,J Z } + ) 

2N(N-1) 
N 2 - 4(J2) _ (J 2 + J 2 ) - N/2 
4N(N- l) ,V 

(Jl) 



v± 



x± 



w 



(6) 



(7) 



N(N- 1) 



N(N — 1)' 



where [A,B]. 

j2 _ NfN 



AB + BA. 



J 2 



Vi for E a=x ,y,z 
2 "(y + !)■ F° r the symmetric states with parity 
conservation, we find x± = from Refs. [O, H2|. Hence 
the pairwise reduced density matrix is shown in X form. 

From the definition in Eq. Q , the quantum discord T> 
and the classical correlation C can be obtained over the 
reduced sub-system von Neumann entropy Ha, the joint 
entropy Ha,b, & n d the conditional entropy H^-^b, which 
are analyzed in detail in Appendix A. 



III. QUANTUM CORRELATION OF THE 
DICKE MODEL 

We study the QD in the Dicke model. The Dicke 
Hamiltonian can be written in terms of the collective mo- 
mentum form 33, 34] 



ickc 



toa) a + AJ Z 



2A . + 



(8) 



where and a are the bosonic annihilation and creation 
operators of the single-mode cavity, A and w are the 
transition frequency of the qubit and the frequency of 
the single bosonic mode, A is the coupling constant. J x 
and J z are the collective spin operators. It is well known 
that this model undergoes a second order quantum phase 
transition from the normal phase to the super-radiant 
phase, separated by the critical point A c = VljK/2. 

We first apply the Holstein-Primakoff transformation 
to change the collective ang ular oper ators to the boso n 
operators b(tf) by J + = bWN - tfb, J_ = V 'N - tfbb, 
and J z = tfb - N/2, where [6,6+] = l [33]. Then the 
displacements of the boson operators are introduced to 
depict the behaviors of super-radiation phase as c+ = a+ + 



^/~Na and & = 6+ — ^/Nj3. By using large N expansions 
of -ffDicko with respect to the new boson operators c+ and 
S up to the I/AT, we obtain the ground state energy as 



E G (a,P) 
N 



U!Ct 



4Aa/Vl- P 2 + A(/3 2 -1). (9) 
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Minimizing the ground state energy gives 



Loa - 2X/3y/l - /3 2 = 
2aA v /l^3 2 "-^£l-/3A = 0. 



(10) 



then we have 



P 2 = max{0,-(l-A 2 /A 2 )}, 
a = ^p^T^. 



(11) 



where the transition point A c = VwA/2. Next we can de- 
rive the matrix elements of the pairwise reduced density 
in Eq. © up to O(l) 

v+ = P\ = (l-/3 2 ) 2 , (12) 
w = y = p 2 {l-p 2 ), u = p 2 {l-p 2 ). 

The von Neumann entropy of the subsystem A and B are 
thus given by 

H(A(B)) = -/3 2 log/3 2 -(l-/3 2 )log(l-/3 2 ), (13) 
H(A,B) = -[2/3 2 (l-/3 2 )]log[2/3 2 (l-/3 2 )] (14) 
-[(/3 4 + (l-/3 2 ) 2 )]log[/3 4 + (l-/3 2 ) 2 ]. 

Note that the quantum mutual information is obtained 
by 1(pa,b) = H(A) + H(B) - H(A,B), so we can derive 
the Eq. Ell and Eq. Ell as 

{X K+ +X k ^) = rj fc (0), (15) 

(X fc;+ -X fe> _) = (2/3 2 - 1)^(0), 

\Y k \ 2 = \Y\ 2 = /? 4 (1 - /3 2 ) 2 sin 2 20cos 2 0, 

Pfc = %(0), 

with fc = 1, 2. 



Therefore the eigenvalues of the conditional density ma- 
trix in Eq. IA3I are obtained as 



2p k 



Vk(0) 
2pk 



{l±[(2/3 2 - V 2 



11/2 



}■ 



By denning x k (0,0) = {(2f3 2 - l) 2 + fj] 1/2 (0<a*<l), 
we simplify the A^(6», 0) = V0±[l±x k (9,<p)]. The condi- 



tional von Neumann entropy is described as 
H(A\{U B })(9^) 
= £{%(0)]og»ft(0) 



fe=l,2 

Me) 

1 2 



(l+x fc (0^))]log[^(l+x fe (0,0))] 



(i-* fc (M))]iog[^(i-a*(M))]} 



log2- ^ 



_ -Kl+z fc (M))log(l+x fe (0,0)) 

fc=l,2 ^ 

+(l-x k (9,<P))log(l-x k (d,cl>))l (17) 

Since = -[(1 + x) log(l + x) + (1 - x) log(l - x)] is 
the monotonically decreasing function, (/) is selected to 
to maximize x k {9,<j))- Furthermore, we find 9 = tt/4 to 
minimize the conditional von Neumann entropy, with 

M = x k (ir/4, 0) = [(2/3 2 - l) 2 + 16/3 4 (1 - /3 2 ) 2 ] 1 / 2 . (18) 

As a result, the conditional von Neumann entropy is 
shown as 

1, 



H(A\Il a ) = log(2) - - [(1 + M) log(l + M) 
+(1 - M) log(l - AO]. 



(19) 



Combining Eqs. CE]), dT4j) , and ([19]), we derive the 
quantum discord as 



V 



-/3 2 log/3 2 -(l-/3 2 )log(l-/3 2 ) 
-[2/3 2 (l - /3 2 )] log[2/3 2 (l - /3 2 )] 
-[/3 4 + (1 - /3 2 ) 2 ] log[/3 4 + (1 - /3 2 ) 2 ] 



(20) 



log 2 

M) + (1 - A/)log(l - M)]. 
While the classical correlation is 



|[(l + M)log(l 



C = -/? 2 log/3 2 -(l-/3 2 )log(l-/3 2 )-log2 (21) 
+ M) log(l + M) + (1 - M) log(l - M)]. 

We next investigate the quantum and classical correla- 
tion in the finite size Dicke model. Two of the present au- 
thors and collaborators have pro pos ed a numerically ex- 
act technique to the Dicke model [34| up to very huge size 
by using the basis of extended coherent states. This effec- 
tive approach has been confirmed recently by comparing 
with the results in terms of basis of the Fock states [35| ■ 
It was demonstrated [35[ that it is very difficult to ob- 
tain convergent results for large number of atoms based 
on usual basis of the Fock states [Hj]. 

In the numerically exact approach [34j |. the wave 
function can be expressed in terms of the basis 
{\<Pn)b ® li) n )} where |j, n) is the Dicke state with 
j = N/2 and \<p n )b ^ s the bosonic extended coherent state 



k=0 



±=(ai+g n ) k e-^-A/ 2 \0) a , (22) 
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FIG. 1: (Color online) (a) Quantum discord, (b) classical 
correlation, and (c) concurrence as a function of the coupling 
constant in the Dicke model with different sizes for A = f.O 
and uj = 1.0. Inset in (a) is the finite size scaling of the 
quantum discord at the critical point A c , and the solid line 
scales as N~ 2/3 . 



where g n = 2An/wV N , N tr is the truncated bosonic 
number in the space of the new operator A n = a + g n , 
\0) a is the vacuum as a\0) a = 0, and the coefficient c n> k 
can be determined through the exact Lanczos diagonal- 
ization. Then, we can derive the elements of pairwisc 
density matrix in Eq. (|6]). Without loss of generality, we 
mainly focus on the resonant case A = uj in the following. 

As is well known that the second-order quantum phase 
transition occurs in the thermodynamic limit. The quan- 
tum correlation is deeply related with the quantum phase 
transition. The reminiscent nature of the quantum phase 
transition in the large size system should be very interest- 
ing. Therefore, we study the size dependence of the QD 
and classical correlation in the Dicke model. In Fig. Q] 
(a) and (b), we displayed the quantum correlation and 
classical correlation as a function of the atom-cavity cou- 
pling constant for several system sizes. The results in the 
thermodynamic limit are also listed. 

It is shown that both the quantum correlation and 
classical correlation are quite small in the normal phase 
(A < A c ). They tend to zero with the increase of the 
atomic number, agreeing well with the results in ther- 
modynamic limit, which is exactly zero. In the ther- 



modynamic limit, there is no excitation of the system 
in the normal phase, so correlation between two arbi- 
trary atoms, which can only be mediated by the pho- 
tons, should be absent. While in the super-radiant phase 
(A > A c ), the QD shows non-monotonous behavior as 
the coupling increases in Fig. Q] (a). The maximum of 
the QD is away from the critical point, even in the ther- 
modynamic limit, locating at 1.28A C . It is different from 
one dimensional XXZ model [l6j , where the maximum of 
the QD is right at the critical point. 

It is very interesting to observe a power law scaling 
T>\\ = \ c ocN^ fl for different detunings, as shown in the in- 
set of Fig. [T] (a). The asymptotic slop in the log- log scale 
in the large size atomic systems for different detunings 
suggests the same exponent fj, = 2/3. To the best of our 
knowledge, such a finite size scaling for the QD itself has 
never been reported in other critical systems. 

As the coupling strength increase further, the QD be- 
comes smaller. In the strong coupling limit, the QD fi- 
nally decreases to zero, which means the non-local quan- 
tum correlation of the system disappears. It can be also 
seen in the Eq. (J2DH that /3 2 = 1/2 as A^oo, so V->0. 

The classical correlation in the super-radiant phase in- 
creases monotonously with the coupling strength, similar 
to those in the different spin chains [lg, l2l| . In the ther- 
modynamical limit, 1 — >l/2 in the deep coupling regime, 
so that H{A) = H(B) = log 2, and H(A,B) = log 2. For 
the conditional von Neumann entropy, H(A\H B ) = 0. 
Hence, the classical correlation saturates at log 2 in the 
strong coupling limit. 

It is well known that the concurrence and the QD are 
both good measures to investigate the nonclassical cor- 
relations. Hence, comparisons between these two quan- 
tities are highly desirable. The scaled concurrence has 
been calculated previously (3^- The concurrence in the 
very large system size has been calculated later by two 
present authors and collaborators [Uj]. We also collect 
the concurrence as a function of coupling constant for 
different size in Fig. [T] (c) for convenience. 

Interestingly, we find that the quantum discord and 
the concurrence show essentially different behaviors in 
the whole coupling regime. First, the concurrence in- 
creases with the atomic number in the normal phase, 
it is always finite in the thermodynamic limit, in sharp 
contrast with the QD in the Dicke model. It is shown 
that the QD decreases with the atomic number in the 
normal phase and vanish exactly in the thermodynamic 
limit. Therefore we conclude that the concurrence is not 
sufficient to quantify the non-classical correlation. Sec- 
ond, the concurrence reaches its maximum at the crit- 
ical point where the QD is still very small. Third, in 
the supperradiant phase, the QD show the broad maxi- 
mum and decay slowly, while the concurrence decreases 
monotonously and rapidly with the coupling constant in 
the strong coupling regime. E.g. for A/A c = 2.0, the 
concurrence becomes negligible, and the value of QD still 
remains about half of the maximum, demonstrating the 
robustness of the QD. In this sense, one can say that 
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FIG. 2: (Color online) First derivative of the quantum dis- 
cord (a) and the classical counterpart (b) as a function of the 
coupling constant in the Dicke model with different system 
sizes for A = 1.0 and uj — 1.0. Inset in (a) shows the finite 
size scaling of the maximum of the dD/dX. 



IV. QUANTUM CORRELATION OF THE LMG 
MODEL 



We turn to the other quantum collective spin model, 
LMG model. Its Hamiltonian reads [26|, [29|, [30] 

ffLMG = - ^ »>; + ^> V i) ~ \ E ( 23 ) 

i<j i 

where Oiii — x,y,z) are Pauli spin-1/2 operators, A is 
magnetic field, 7 is anisotropic parameter. In the frame- 
work of the collective spin operators = J^i °f 7 2 with 
k = x, y, z, and J± = J x ±iJ y , the model can be rewritten 
as 

H LMG = -AJ,-l[J 2 + 7 J 2 -/V(l + 7 )/4].(24) 

We focus on the case of 0< 7 < 1 and A > where the 
second-order quantum phase transition can occurs. 

We also first study the quantum and classical cor- 
relation in the thermodynamic limit. Based on the 
Holstein-Primakoff Transformation J z = b^b — N/2, J + = 
tfy/N — Wb, J_ = y/N — tfbb, We displace the boson op- 
erator to the form of = b^ + a. By large TV expansion 
of i?LMG ; the ground state energy per spin is shown as 

E G /N = -[(1 -a 2 )a 2 + \{a 2 - 1/2)]. (25) 

By minimizing Eg, we find 



QD is a pretty good measure of the quantum correlation 
independent of entanglement. 



To quantify the quantum phase transition of the Dicke 
model, the first derivative of the QD and the classical 
correlation as a function of the coupling constant are 
presented in Fig. [2] The cusp-like peak of 4? curves 
emerges with increase of the atomic number, indicating 
an non-analytical behavior. The peak position tends to 
the critical point with increasing N. It would provide 
a good tool to detect the critical point by a finite size 
study. In the spin chains, the first derivative of the QD 
also shows its peak right at the critical point, which is 
also non-analytical and discontinuous. 



Interestingly, the logarithmic scaling of the maximum 
of the dD/dX is also observed in the large size regime, 
which can be fitted well by (^) max = 0.377 log 2 (JV) - 
0.401, demonstrated in the inset of Fig. [5J This loga- 
rithmic divergence has been also reported in Ref. [l7], H|[ 
for XY spin chain and transverse field Ising chain. The 
first derivative of the classical correlation ^ display a 
different behavior. It becomes more rounded around the 
maximum with the increase of the size. 



a 2 = min{l, 



1 + A 



}■ 



(26) 



Consequently, we derive (J z )/N = (a 2 -1/2), (J 2 )/N 2 = 
(a 2 - 1/2) 2 , and {J\)/N 2 = a 2 (l - a 2 ). Then the el- 
ements of the pairwise density matrix in Eq. ([6]) up to 
O(l) are derived by 



-) 2 , »- = ( 



1 - A, 



1 + A 2 

2 ' ' " v 2 



(27) 



1+A 1-A I 

w = y = {—^— )(— )> u = {— 



-v^>. 



Similar to the expressions of the elements in Dicke 
model in Eq. (fl"2"j). we obtain the corresponding relation 
in the LMG model 



1 
(3 2 



1-A 
1 + A 



(28) 



then the elements can expressed as 

v + = (1-/3 2 ) 2 , 
v- = P\ 

□2 



(29) 



w 
u 



y = ^(l-/3 2 ), 
/3 2 (l-/3 2 ). 
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FIG. 4: (Color online) First order derivative of the quantum 
discord (a) and the classical counterpart (b) as a function of 
the field in the LMG model with different system sizes for 
7 = 0. Inset in (a) exhibits the finite size scaling of the 
minimum of the dT)/d\. 



FIG. 3: (Color online) (a) Quantum discord, (b) classical 
correlation, and (c) concurrence as a function of the field in 
the LMG model with different system sizes for 7 = 0. Inset 
in (a) is the finite size scaling of the quantum discord at the 
critical point A c , and the solid line scales as N 2 ^ 3 . 

The influence of the anisotropic parameter 7 on the 
quantum correlation is checked, and we find that the QD 
is almost not altered with 7. Hence we mainly focus on 
7 = in the following. 

The finite size LMG model can be studied by the exact 
diagonalization on the basis of the collective spin opera- 
tors [2{|. For the ground-state, the convergent results can 
be obtained for very huge system. The size dependence 
of the QD and classical correlation is given in Fig. [3] The 
QD shows the non-monotonous behavior in the symme- 
try broken phase (A < A c ). Starting at 0, the QD then 
shows the single broad maximum at A = 0.62A C , far from 
the critical point. While the QD is very small in the sym- 
metry phase (A > A c ). The classical correlation decreases 
monotonously in the symmetry broken phase, approach- 
ing to zero at the critical point. If the magnetic field 
A = 0, the classical correlation is easily obtained as log 2. 
While in the symmetry phase, the classical correlation is 
rather small for finite size system and becomes zero in 
the thermodynamic limit. 

Similar to the Dicke model, the QD of the LMG 
model also shows the power law scaling obviously as 
T>\ = \ c (xN~ fJ - at the critical point, as shown in the in- 
set of Fig. |3fa). The scaling exponent fi in the large N 
regime is very close to 2/3, the same as that obtained 



above in the Dicke model, providing a new piece of evi- 
dence of the same universality class of these two models. 

The QD in the LMG model in the thermodynamic limit 
has been preliminarily studied in Ref. [17| . By mapping 
LMG model to the two bands fermion model exactly, 
they alternatively studied the quantum correlation and 
the classical counterpart of the fermions based on the 
corresponding density matrix. They found that both the 
QD and the classical counterpart decrease monotonously 
in the symmetry broken phase with the magnetic field, 
and keep zero in the symmetry broken phase. In the 
present paper, the pairwise density matrix is directly de- 
rived from the qubit in LMG model, which is related 
but not identical with that in Ref. [17| • The different 
definition may account for the different QD obtained in 
two papers. Our results for the classical correlation are 
consistent qualitatively with those in Ref. [17] . 

The comparisons of the QD and the scaled concurrence 
in the LMG model are also made. The scaled concur- 
rence for large system size has been calculated by Dusucl 
et al[37j], which was shown in their Fig. 9. For complete- 
ness, we recalculate the concurrence for much too larger 
system size and list them in Fig. OUJc). At the critical 
point, the concurrence shows the maximum, while the 
QD becomes 0. In the symmetry phase (A > A c ), the 
quantum discord are very small, implying no quantum 
correlation. While the concurrence remain finite and de- 
crease gradually and monotonously with the field A. 

In order to study the critical behaviors of the LMG 
model, we display the first derivative of the QD and clas- 
sical correlation in Fig. 4? has a pronounced mini- 
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mum around the critical point. Note that the minimum 
position approaches the critical point, providing a con- 
vincing method to locate the critical point. This behavior 
strongly implies the sudden transition from the symme- 
try broken phase to the symmetry phase. The finite size 
scaling of the minimum of the dV/dX is also performed in 
inset of the Fig. [3] The logarithmic divergence fitted by 
^| min = -0.044 log 2 (iV) - 0.346 is also well exhibited in 
the large size regime, similar to that in the Dicke model. 
It should be point out that the coefficients of log 2 (N) in 
the two models are different. Combined with the similar 
scaling behavior reported in spin chain [13, HH, such a 
logarithmic divergence of the first derivative of the QD 
may be universal in critical systems, further confirma- 
tions in extensive studies are however needed. The clas- 
sical correlation shows the broad valley around 0.55A C , 
and becomes as the field parameter A > A c . 



V. SUMMARY 

The quantum discord and its classical counterpart in 
the Dicke model and the LMG model have been inves- 
tigated in this paper. These correlations in the ther- 
modynamic limit have been derived analytically and in 
the finite size system are obtained numerically. Power 
law scaling behavior at the critical point is observed in 
both models. Such a scaling behavior has not been re- 
ported in other critical systems, as far as we know. The 
same scaling exponents provide a new evidence of the 
same universality of two models. The position of the 
pronounced maximum or minimum of the first derivative 
of the QD approaches to the critical point with the in- 
crease of the system size, which can be used to identify 
the quantum phase transitions. A logarithmic diverging 
behavior for the first derivative of the quantum discord 
is also observed in the two models. The coefficient of the 
logarithmic term is model dependent. 

We also find that the quantum discord and the concur- 
rence show essentially different behaviors in both models. 
In the normal phase of the Dicke model and the sym- 
metry phase of the LMG model, the QD is zero in the 
thermodynamical limit and extremely small in the finite 
size system, while the concurrence remain finite. In the 
supperradiant phase of the Dicke model and the symme- 
try broken phase of the LMG model, the QD show the 
broad maximum and decay slowly, while the concurrence 
decreases monotonously and rapidly when away from the 
critical point. These salient features might be used for 
the applications of the quantum correlation in some ex- 
perimentally realized systems to the quantum informa- 
tion science and the quantum computing. 
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Appendix A: Derivation of the quantum correlation 

Considering the pairwise atom reduced density matrix 
in Eq. |B1 the reduced density matrix p A is obtained by 

PA = Ttb{pa,b} 

= {v++w)\e) A {e\ + (v-+w)\g) A {g\. 

The Von Neumann entropy of the reduced system A is 
shown as 

H{A) = — (v + + w) log(w+ + to) — (v- + w) log(ti_ + to). 

(Al) 

The joint von Neumann entropy can be derived from 
the joint pairwise matrix in Eq. [5] By solving the Eq. [51 
we obtained 

(A + y - w)[(w - A) + y}[(\ - v+)(X - v.) - \u\ 2 \ = 0. 

The eigenvalues — 1, 2, 3, 4) can be obtained analyt- 
ically. Considering x± = 0, the joint entropy is demon- 
strated as 

H(A,B) = -(to + y)log(to +y) - (to - y)log(u> - y) 

- Alo § A ' ( A2 ) 



A=A± 



where 



A± = + v-)±[(v+ - v-) 2 + 4M 2 ] 1/2 }- 

The conditional density PA\n® i s measured by the pro- 
jections tuned by 9 and 4> 

|*!) B = cos(0)|e) B W^ S m(0)|<7> B 
|* a ) B = e -^sin(0)| e ) B -cos(0)|<7) B . 

Under such projections, the conditional density matrix is 
shown as 

PA\n% a = l")eH{l e M e l X ",+ + \9)a{9\ x c- 

+\e) A (g\Y a + \g) A (e\Y*}/p a . (A3) 

For a = 1, 

X h+ = v+cos 2 (6) + wsm 2 (6), (A4) 
= wcos 2 (6>) +i;_sin 2 (6»), 

pi = to + v + cos 2 (8) + vsm 2 (9). 
For a = 2, 

X 2 ,+ = v + sm 2 (8)+wcos 2 (8), (A5) 
X 2 ,- = wsin 2 {8) +t;_cos 2 (6»), 

Y 2 = -sin(60cos(6>)[e I< V + e -i *y], 
p 2 = W + v + sm 2 (9) + v-Cos 2 (9). 
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Then the eigenvalues of the density matrix Eq. IA3I read 

"Pot 

+4|F a | 2 ] 1/2 }- (A6) 

Finally, the conditional Von Neumann entropy is shown 
as 

H(A\{n®}) = - P«[W(«logA£(M)(A7) 

01=1,2 

+X1(e,0ogX1(e, </>)}. 

Hence 

5(0,0) = H(A)-H(A,B) + H(A\Tl B ). (A8) 



The quantum discord can be obtained by optimizing the 
9 and <j> both in the regime [0,7r/2] to minimize 6(9, 4>), 
shown as 

V = min {6(9, (p)}. (A9) 

Consequently, the corresponding classical correlation can 
also be obtained by 

C = nmx{H(A) + H(B) - H(A, B) - 6(9, </>)}. (A10) 

{8,4>} 
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